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Abstract. In this paper we have examined the linear stability of triangular equilibrium points in the 
generalised photogravitational restricted three body problem with Poynting-Robertson drag. We have found 
the position of triangular equilibrium points of our problem. The problem is generalised in the sense that 
smaller primary is supposed to be an oblate spheroid. The bigger primary is considered as radiating. The 
equations of motion are affected by radiation pressure force, oblateness and P-R drag. All classical results 
involving photogravitational and oblateness in restricted three body problem may be verified from this 
result. With the help of characteristic equation, we discussed the stability. Finally we conclude that 
triangular equilibrium points are unstable. 
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1. Introduction 

We study the motion of three finite bodies in the three body problem. The problem is restricted in the 
sense that one of the three masses is taken to be so small that the gravitational effect on the other masses 
by third mass is negligible. The smaller body is known as infinitesimal mass and remaining two massive 
bodies as finite masses. The classical restricted three body problem is generalised to include the force of 
radiation pressure, the Poynting-Robertson effect and oblateness effect. J.H. Poynting(f 903) considered the 
effect of the absorption and subsequent re-emission of sunlight by small isolated particles in the solar system. 
His work was later modified by H.P. Robertson(f 937) who used a precise relativistic treatments of the first 
order in the ratio of the velocity of the particle to that of light. The motion of an infinitesimal body near 
one of the Lagrangian points is said to be stable if given a very small displacement and small velocity, the 
body oscillates for a considerable time around the point. If it departs from the point as the time increases, 
the motion is said to be unstable. Mathematically, we express that the stability of motion means bounded 
displacement and velocity with bounded functions of time in the neighborhood of the equilibrium points. 

If a satellite flies high enough above the Earth and is large enough in size, but at the same time has a 
sufficiently small mass, then the pressure of the Sun's light rays (known as radiation pressure ) has a very 
strong effect on its motion. The distance of the satellite to the Sun practically does not change, and so the 
magnitude of the radiation pressure force is practically constant. The Poynting-Robertson (P-R) effect is 
the most important nongravitational and nondisipative interplanetary dust particles. The P-R effect may 
be significant in dealing with other problems in astrophysics. 
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The effect of radiation pressure and P-R drag in the restricted three body problem has been studied 
by Colombo et aZ.(1966),Chernikov Yu.A.(1970) and Schuerman(1980) discussed the position as well as the 
stability of the Lagrangian equilibrium points when radiation pressure, P-R drag force are included. Murray 
CD. (1994) systematically discussed the dynamical effect of general drag in the planar circular restricted 
three body problem. Liou J.C. et aZ.(1995) examined the effect of radiation pressure, P-R drag and solar 
wind drag in the restricted three body problem. The classical three body problem has five Lagrangian points. 
Their location and stability properties are well known. The three collincar points Lx^L^Lz lies on the line 
joining the two primaries, are linearly unstable. The two Lagrangian points L4, L§ lie in the orbit of the test 
particle which are symmetrical about the line joining the primaries, are linearly stable with the condition 
that the mass parameter /j, to be less than .0385. Hence, we are interested to discuss the linear stability of 
triangular equilibrium points in generalised photogravitational restricted three body problem with Poynting- 
Robertson drag. After Newton solved the problem of the orbit of single planet around the Sun, the natural 
next challenge was to find the solution of two planets orbiting the Sun(i.e. three body problem ). Many of 
the best minds in Mathematics and Physics worked on this problem in the last century. Three-Body problem 
is a continuous source of study, since the discovery of its non-integrability due to Poincare(1892). Regular 
and chaotic motion have been widely investigated with any kind of tools, from analytical result to numerical 
explorations. The restricted three-body problem is a modified three-body problem when one of the three 
masses is set to zero(infinitesimal mass ). 

The restricted three body model, Earth-Moon and Spacecraft is insoluble, although specific solutions 
exist, like the ones in which the spacecraft is positioned at one of the Lagrangian points. There are many 
periodic orbits in the restricted three body problem. One of the most famous, discovered by Lagrange, is 
formed by an equilateral triangle. In "Earth-Moon-Space Station model", if the Moon is not too massive, 
this orbit is thought to be stable. If the space station is pushed a bit to one side (in position or velocity ), 
it's supposed to make small oscillations around this orbit. The two kinds of triangular points are called L4 
and L5 points. 

There is a group of enthusiasts who want to setup a colony at L5 point of the Earth-Moon system. As 
already noted, because L4 and -L5 are the stable points of equilibrium, they have been proposed for sites of 
large self-contained "Space colonies", an idea developed and advocated by the late Geraid O' Neill (1974). 
The three body problem have an interesting application for artificial satellites and future space colonization. 
Triangular points of the Sun- Jupiter or Sun-Earth system would be convenient sites to locate future space 
colonics. Such colonies would presumably mine asteroids for at least some of their raw materials. Present day 
technology permits the construction of spacecraft whose value of /3(ratio of radiation force and gravitation 
force) can be made or varied in the range < (i < 1. (A spacecraft with /3 > 1 has been considered by 
NASA as a "Sailing" vehicle for chasing Hallcy's comet [Friedman 1976]). A Space Colony with a solar fac- 
ing adjustable "sail" could thus "park" (i.e. remain stationary with respect to the Sun and a planet without 
expending large amounts of energy) not only at the classical positions but at any heleocentricdistance out 
to that of the planet. 

Thus we have included the purpose of the work and its significance in the field of Celestial Mechanics 
and Dynamical Astronomy. Application of results to realistic actual problem is obvious. L4 and L 5 points 
became unstable due to P-R drag which is very remarkable and important, where as they are stable in 
classical result. 
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2. Location of Triangular Equilibrium Points 
Equations of motion are, 

(1) x-2ny = U x , where, U x = ^--—^- 

ax rf 

(2) y + 2nx = U y , U y = ^ - 

(3) z — U z, U z 



dz r{ 



(4) m = 



n 2 (x 2 +y 2 ) (1-m)«i , A* , M 



2 



ri r 2 2r| 



r 2 = (x + M ) 2 + y 2 , r2 = (x + M -l) 2 + y 2 , n 2 = 1 + ^A 2 , 

(x + li)N . yN . , zN 

Ni = 5 V x — ny, A 2 = '-^r + y + n(x + u), N3 — —5- + z, 

r{ rf ' rf 

N = (x + fi)x + yy + zz 

r 2 — r 2 

^4-2 = " 5r ^ p be the oblateness coefficient, r e andr p be the equatorial and polar radii respectively r be the 
distance between primaries, q = (1 - J?) be the mass reduction factor expressed in terms of the particle's 

radius a, density p and radiation pressure efficiency factor \ ( m the C.G.S. system) i.e., g = 1 — 5,6 — 
Assumption g = constant is equivalent to neglecting fluctuation in the beam of solar radiation and the effect 
of solar radiation, the effect of the planet's shadow, obviously q < 1. Triangular equilibrium points are given 
by U x = 0, U y — 0, z — 0, y 7^ 0, then we have 



f nWi[(l-/i)(l + |A 2 )+M(l-#) 
(5) = xoll v " ' 



s 2 




(6) y* = yo< 1 



3/x(l - p)y x 

^<5 2 [2/x - 1 - /i(l - + 7(1 - ^ 2 (1 - 4^)A 2 



1/2 



where xq = 4 — //,yo = ±<5(l — ^-) 1 ^ 2 and 5 = g^ 3 



3//(l - p)yg yg j 

1/3 
= ?i 

where subscript * indicate the value at the equilibrium points. Eq.(|3J and (J^J are valid for W\ <C 1, 
A 2 < 1. 

3. Stability of Triangular Equilibrium Points 

We suppose a = Aie xt ,(3 = Bie xt ,-f = C\e xt be the small displacements from Lagrangian points (x», y», z») 
and Ai, Bi,C\, A are parameters, x = x* + a, y — y» + (3 and z = z* + 7. Therefore the equations of 
perturbed motion corresponding to the system of Eq.(l),(2),(3) may be written as, 

(7) a - 2n/3 = U° x + aU° xx + pU° xy + 7 [& + aU° xx + pU° xy + W xi 

(8) (3 + 2nd = U a y + aU° yx + pU° y + 7 E7j, + aU° yx + 0U° y + jUyi 

(9) 7 = U° z + aU° zx + (3U° zy + 7 U° ZZ + aU° zx + (iU% + 

At (s* U° x = C/° = C/° = l& = U°, = U° zx = U% = U° xz = U yz = U° zx = U° z = 

Above system of equations can be written as, 

(10) (A 2 - AE& - U° xx ) + [-(2,i + U° xy )\ - V%]Bi = 

(11) [(2n - £^)A - [/, J^i + (A 2 - \U° yy - U° yy )B x = 

(12) (A 2 - \U° ZZ - f/°Jd = 
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Eq.(H3 gives us, A = u °^ ± \ / ^ 2 + 4 Hl±_ From Eq .( 3 ) 



u° zz = -f, u° z 



Wi 



where, / 



(1 - n)qi \i 



3^2 

2r| 



(13) 



A = 



lj Wi ± /Wf 



4/ } as Wi > and / > 



Since A consists always negative real part, so motion is asymptotically stable in z direction. Now Ea. l|10|) 
and f lllfl has singular solution if, 



X 2 — xu2a, — u2. 



XX XX 



(2n - UL)X - U° x 



-(2n + U^)X-U% 
X2 ~ XU yy ~ U yy 



= 



(14) => A 4 + a\ 3 + b\ 2 + cX + d = 

At the equilibrium points the system of Eq.(l),(2) gives 



(15) a = 3-o— , o = 2n — /* 



^5 



6 + &i &o = 2n ,6i = -/* 



3/xA 2 2Wf 



(16) 
(17) 

where 



d = (« 2 - /.) 



M2 

n 2 + 2/* - 



2 r 5 
1*'2 

3juA 2 



9/i(l - fi)qi / 5^2 \ 2 
V 2r3 / y * 



2^ 



6nWl 



1*' 2i 



2* • 



(1 + ^r-) {(x* + - (z* + ft) + vl) + 
nWi 



n 2 W 2 
"1* 



!/3 j -1 



M\ ^ _ 1 | nWi 

6(1 - M)yo 2 J ' 2 * 3/ij/o 



1 - -A 2 
2 



Let (7 = l — e with e > very small real quantity then 

n^ 1 (9 + 2e)(4-3A 2 ) (1-^2 
'* 18V3 + 2 

We can write the four roots of the classical characteristic equation as, 

(18) X j =±zi, where, z 2 = i{l T [l - 27 M (1 - m)] V2 }, J = 1,2,3,4. 

Due to P-R drag and oblateness,we assume the solution of the characteristic Eq. ( I14f> of the form 

A = Aj(l + ei + e 2 i) = ±[— e 2 + (1 + ei)i]z 
where ei, e 2 are small real quantities. We have in first order approximation 



A 2 = [-(1 + 2ei) — 2e 2 i]z 2 , A 3 = ±[3e 2 - (1 + 3ei)«]z 3 , A 4 = [(1 + 4ei) + 4e 



Substituting these values in characteristic Ea.f ll4[) neglecting products of e\ or e 2 with a, 61, c, we get 

— d + b a z 2 — z 4 



(19) 
(20) 



ei 



a 2 



2z 2 (2z 2 - bo) 
=fcz ± az 3 
2z 2 (2z 2 - 6 ) 



If e2 7^ 0,then the resulting motion of particle displaced from equilibrium points is asymptotically stable 
only when all the real parts of A are negative. For stability, we require Re(X) < 0. Rate of growth or decay 
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is determined by the largest real part of the roots of the characteristic equation, where Re(X) denotes the 
real parts of A and given by 



(21) 

we have, 
(22) 

Taking positive sign 



Re{x) = W^T) 

i{l T [l-27Ml-M)] 1/2 } 
27 

z 2 = 1 - - M (l - fl) 



Consider 



r 27 if 27 i _1 

i?e(A)<0 then |c- a[l - — M (l - /x)] )[l - — M (l - M )j < 



(23) 

Taking negative sign, 
(24) 



27 

c + — A 4 ) (2c — a) < a 



r27 i r 27 

Re(X) = - |c - a—fi(l - A*) J [1 - — - /u) 

27 

< c+ — - |i) (2c- a) 



< 



27 

< c + -yAHl _ M)(2c- a) < a 



From above inequalities we get, 

(25) 

(26) 0<c<a as ^ ^ Murray CD. (1994) 

Inequality (ESI is necessary condition for stability of triangular equilibrium points L4 and L5. 

4. Conclusion 

(i) When Wi = 0, then 

a = 0,c=0,b = n 2 -3fiA 2 , 



Equation f 114(1 gives us, 



d=9/i(l-At) 1 + 77^2 1- 



nV2 



2/3 
gl 

4n 4 



x2 -b±(b 2 -Ad) 
X = ^— For stable motion b > Ad 



2/3 
1x \ 2 

— - — n 



i.e., (n 2 -3A) 2 >36//(l-/x)(l + -A 2 )(l ^ ]; . 

(ii) For classical case A 2 = 0, gi = 1, W\ — 0, n = 1. Then from inequality ( 126(1 we have 

1 > 27ai(1 - At) => At < -0385 

(iii) When Wt ^ 0, A 2 = 0, n = 1. Then 

2VF 2 



a = 3^—, & = 2 — /* 



I r 2 rf J 



^(i-/0(i + 2/,) + ^^^-r&[(^+M) 2 -(^ + At) + ^ + ¥ 



r 2 r\ 



as in Chernikove Yu.A.(1970) 
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If we consider the P-R drag force then c < 0. This dose not satisfy the necessary condition of stability. 
Hence motion is unstable in the linear sense. 
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